Wannier- functions characterization 
of floating bonds in a-Si 

M. Fornari/^'^) N. Marzari/^^ M. Peressi*^^) and A. Baldereschi^^'^^ 

Naval Research Laboratory, Washington, DC 20375 and CSI George Mason 

University, Fairfax, VA 22030 
(2) IJ^PM and Dipartimento di Fisica Teorica - Universitd di Trieste, Strada Costiera 

11, I- 3 40 14 Trieste, Italy 
^■^^ Institut de Physique Appliquee, Ecole Polytechnique Federale de Lausanne 
PHB-Ecuhlens, CH-1015 Lausanne, Switzerland 

Abstract 

We investigate the electronic structure of over-coordinated defects in amor- 
phous siUcon via density-functional total-energy calculations, with the aim of 
understanding the relationship between topological and electronic properties on 
a microscopic scale. Maximally-localized Wannier functions are computed in 
order to characterize the bonding and the electronic properties of these defects. 
The five-fold coordination defects give rise to delocalized states extending over 
several nearest neighbors, and therefore to very polarizable bonds and anoma- 
lously high Born effective charges for the defective atoms. 

1 Introduction 

Coordination is a well-defined concept in crystalline systems, where it is closely related 
to the concept of bonding. For instance, silicon atoms are tetrahedrally coordinated, 
and exhibit covalent bonds with their four nearest-neighbors (NNs) in the ordinary 
crystalline phase (c-Si). Whenever symmetry is broken, as in the case of amorphous 
silicon (a-Si) studied here, the link between coordination and bonding becomes less 
obvious. Even if atomic positions are intimately connected with the electronic ground 
state, we are usually unable to extract detailed informations about bonding from 
a purely geometrical investigation. Standard analytical tools based on coordination- 
number analysis are thus insufficient to characterize bonding properties, since they are 
insensitive to the details of the electronic charge distribution. Whenever strong topo- 
logical disorder is present, as in the case of atoms with three or five nearest neighbors 
(T3 or T5 defects), the ionic potential and the charge distribution can be significantly 
different from the crystalline case, and great care should be used to identify the nature 
of bonding. 

In this work we use density-functional theory in the local density approximation 
to calculate the electronic properties of over-coordinated defects in a-Si. Density- 
functional theory provides an accurate and parameter-free description of the electronic 
ground-state, and is equally capable of dealing with four-fold coordinated atoms and 
with more complex topologies. The characterization of defects is then based on the 
decomposition of the electronic ground state into localized orbitals, using the technique 
of maximally-localized Wannier functions (MLWF) 0] . In this approach, the extended 
Bloch orbitals are transformed via unitary transformations into a representation where 



they are maximally localized. The role and importance of such localized Wannier 
functions in the study of disordered systems (amorphous silicon in particular) has been 
advocated by Silvestrelli and coworkers who have shown that the coordination 
analysis is often insensitive to the electronic charge distribution and that similarly 
coordinated atoms can be surrounded by rather different bonding environments. Our 
present results support these findings. 

The conjecture that over-coordination could play an important role in the formation 
of mid-gap electronic levels in amorphous silicon has been recently validated by some 
of us [§], using accurate density functional calculations. In that work we argued that 
T5 defects can be responsible, as much as T3 ones, for states close to the Fermi level. 
We also showed that in the case of T3 defects the mid-gap electronic state originates 
from the dangling bond and is well localized on the T3 defect itself. Conversely, in 
T5 defects the midgap electronic state is delocalized over several NNs of the T5 site, 
in agreement with the findings of tight-binding calculations reported in Ref. [Q. In 
this work we pursue further the study of the extended nature of T5 defects. Our 
goal is to understand how many atoms can be actively involved in a defect, since this 
number influences the shape of the "super-hyperfine" structure of the D center in the 
electron-spin resonance (EPR) spectrum In turn, this should discriminate dangling 
bonds from floating ones. It should be noted that both T3 and T5 centers generate a 
D signal; it is the contributions of the secondary ions involved in the uncompensated 
spin distribution (the super-hyperfine structure) that should be very different in the 
case of a dangling or a floating bond. 

The plan of the work is as follows: in Sec. 2 we describe our ab-initio calculations 
and we review the method used to compute the MLWFs. We also comment on the link 
with the dynamical Born effective charges Z. In Sec. 3 we discuss the results for our 
a-Si sample, in comparison with the results obtained using the "electron-localization 
function" (ELF) and the "atomic-projected charge" (APC) analysis 0. Sec. 4 is 
devoted to the conclusions. 

2 Computational tools 

The calculations performed here are based on density-functional theory in the local- 
density approximation, using a cubic super-cell containing 64 Si atoms 0, a plane wave 
basis set, and accurate norm-conserving pseudopotentials. The electronic ground-state 
is obtained using an all-bands conjugate gradient minimization p]. The sampling in 
reciprocal space is performed with 8 k points in the full Brillouin zone. 

The analysis of the Kohn-Sham orbitals and their charge density has been previ- 
ously performed using the APC and the ELF analysis ^. APC is a measure of the 
contribution of the different atoms to the charge density, and it is obtained projecting 
the Bloch orbitals "^nk on localized atomic functions [Q. ELF, defined as in Ref. [Q, 
is a measure of the conditional probability of having one electron close to another 
one with the same spin. It approaches its upper limit (1.0) when the electron density 
resembles a covalent bond or in the presence of unpaired electrons. 

Both these approaches provide useful results in the characterization of the elec- 
tronic distribution; still, there are several limitations that will be mentioned in the 
next section. The analysis based on localized Wannier functions overcomes these lim- 



Figure 1: Left panel: structure of the "defect group" embedded in our a-Si sample. Atoms 
are labelled T5 (T4) when they are five-fold (four-fold) coordinated. Atom I (interstitial) 
is the only four- fold atom that is a NN of both Ts's. Right panel: the "defect group" is 
shown together with the centers of the associated Wannier functions. CI and C2 are the 
centers with the lowest symmetry, and exhibit the strongest deviations from the geometrical 
mid-bond positions. 

itations; as a byproduct, it also provides informations on the dielectric properties of 
the system. 

Wannier functions are an alternative representation of the electronic ground state, 
that classifies states using a spatial coordinate R and a band index n. Wannier func- 
tions are not eigenstates of the Kohn-Sham Hamiltonian, but are related to them via 
the unitary transformation 

= 777^ / |vI/„k)e-*^ Vk. (1) 
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In the electronic structure problem there are degrees of freedom that do not affect 
the self-consistent ground state, but determine the shape of the WFs obtained from 
Eq. (|l]). They consist of arbitrary unitary rotations UmnO^) that mix together at any 
given k in the Brillouin zone fully occupied Bloch orbitals (in the case of a single band, 
these unitary matrices reduce to a phase factor 0(k)). It is of paramount importance, 
in order to provide a meaningful real-space representation, to choose these arbitrary 
matrices f/mn(k) so that the resulting WFs are well localized. The approach used here 
follows the lines of Ref . [|T| , where the unitary rotations are refined until the resulting 
Wannier functions are maximally-localized, i.e. they have minimum spread 

n = Y,[{0n\r^\0n) - {On\r\On)^]. (2) 

n 

Incidentally, the sum of the Wannier function centers (WFCs) = (0?T,|r|Or;,) is 
directly related to the macroscopic polarization of the sample ||10[, and this makes 



the Wannier function analysis attractive to study dielectric properties. In particular. 




Figure 2: Wannier functions centered near crystal-like coordinated atoms in the a-Si su- 
percell. In each case a covalent bond (light gray) exists between two four-fold coordinated 
atoms, and the WFC is close to the geometrical bond center. For each atom, clover back- 
bonds (dark gray) point towards the remaining three neighbors. 

the change in polarization induced by the displacement At^ of an atom N is directly 
related to the electronic component of its Born dynamical-charge tensor Zn by 
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(3) 



where n runs over all the occupied WFs in the unit cell and Ar^ are the displacements 
of the WFCs. 



3 Results on a- Si 



The presence of mid-gap states in a-Si is usually ascribed to the existence of T3 defects 
in the system. Such conclusion relies on EPR experiments, which measure the unpaired 
spin oscillations of the electrons in the dangling bonds. Nevertheless, certain features of 
the EPR signals (the "super-hyperfine structure" ) derive also from the ionic magnetic 
moment and the quadrupolar term of the secondary ions in the defect Thus, the 



EPR signal is indirectly sensitive to the presence of complex geometries. T5 defects 
are expected to give rise to more delocalized states involving several atoms 0, 
and could thus provide a distinctive shape to the EPR signal. Prompted by these 
considerations, we investigate in this work the electronic properties of over-coordinated 
defects in amorphous silicon. 

Our samples have been obtained by annealing configurations spanned during a 
first-principles molecular dynamics simulation 0. In particular, we have selected for 
this study a sample presenting two T5 coordination defects. Although this corresponds 
to a density of defects higher than that of experiments, it still results in a useful model 
to investigate the effects of local strain and topological disorder. The defects were 
first identified using a geometrical analysis: an atom is counted as a neighbor if it 
resides inside a sphere of radius Rc = 2.6 A (to compare with the theoretical NN 
distance in c-Si (Inn = 2.357 A). This radius is chosen as the first minimum in the 
correlation function g{r), since the peaks correspond, intuitively, to the succession of 
shells surrounding the reference atom. The two T5 defects (labelled T^ and T^ in 
Fig. P belong to a "defect group" (Fig.|l], left panel), mostly composed of bonds which 
are longer and weaker than four-fold ones. Such environment seems favorable to induce 
delocalization of the wavefunction over several neighbors; our previous APC analysis 
(see Ref. 0) supported this conjecture, as did a similar investigation (Ref. [|]) based 
on the electronic localization function. Both approaches however suffer from inherent 
limitations: the projection technique used for the APC does not provide meaningful 
informations on the mid-bond region; on the other side, ELF (once it has been been 
used to unambiguously distinguish between a T3 and a T5 site) is not able to improve 
the qualitative description of the defect that can be obtained from the direct inspection 
of the charge density. 

To improve our understanding of the electronic structure of the defect group, we 
calculated the MLWFs for this configuration. In crystalline Si each WF is centered in 
the middle of the bond and oriented along the bond ; a covalent "barrel" of charge is 
shared between the two bonded atoms, while some back-bonding charge is distributed 
between each of these atoms and its remaining three NNs. We expect that in a slightly 
distorted environment, where short range order is conserved, the shape of the WFs 
would not drastically change. This is clearly confirmed by our calculations (see Fig. 

Small deviations from the crystalline case are due to the bending and stretching 
of bonds, but no doubt exists on the nature of the bond, that resembles closely to the 
bond of crystalline silicon. Such "regular" bonds are clearly related to their crystalline 
counterparts by having a very peaked distribution in their spreads, with a maximum 
around 2.30 (vs. 2.04 for crystalline silicon, using the same 8 k-point sampling). 

We started our analysis focusing on the position of the WFC, as suggested by Sil- 
vestrelli et al. . In that work it was observed that the centers of charge for anomalous 
WFs tend to be closer to one specific atom than those of regular WFs, which are in- 
stead equally distant from two bonded atoms. In Fig.[I| we show the WFCs belonging 
to our "defect group"; it can be seen that around the interstitial atom I the Wannier 
centers depart strongly from their ideal mid-bond positions. Similarly, the shape of the 
corresponding Wannier functions is clearly anomalous, with a delocalization extending 
over more than two atoms and with a much larger spread than average, of the order 
of 4 to 8 A^. This confirms our conjecture of delocalized orbitals connecting different 





Figure 3: Wannier functions corresponding to two bonds included in the "defect group" of 
Fig.|l| we show the WF centered at CI (left panel) and C2 (right panel). CI occurs between 
two bonds, and the corresponding WF is delocalized over the T5 site and two of its nearest 
neighbors. 

atoms inside the "defect group" ; an inspection of the shape of the WF clearly confirms 
the delocalization of the electronic states. In Fig.|^ the left panel shows the WF with 
the center close to the Tf defect (CI); its shape suggest a bond shared between two 
NN. In the right panel of Fig.|^ the delocalized orbital centered at C2 is shown. The 
shape of this orbital is also influenced by the additional interaction between the two 
T5 defects. 

The anomalous shape and spread of some WFs present in our sample suggest that 
unusual polarization properties could be present, and could be related to topological 
disorder. For this reason, we investigated the dynamical charge tensor Z, to quantify 
the local anisotropy and deviation from the crystalline order. Our preliminary results 
suggest that disorder plays a considerable effect on the dynamical charges of a-Si. 

0, 



Following Ref. 



we decompose the Z tensor into an isotropic contribution 
(corresponding to the / = spatial rotation representation), a / = 1 antisymmetric 
one , and a / = 2 traceless symmetric contribution. In crystalline Si the effective 
charges are zero, and the electronic contribution (-4 times the identity matrix) cancels 
exactly the ionic contribution (+4 times the identity matrix). This is not the case in 
amorphous silicon, and we find that atoms belonging to topological defects exhibit a 
very different behavior. To this purpose it is instructive to inspect Table |I|, where the 
electronic effective charge tensors of some selected atoms are given, together with their 



decompositions. Even for regular T4 atoms (T^ and 



the effective charges can 



show strong anomalies, although this could be an artifact of the high density of defects 
in the sample. The anisotropy around atom T^ on the other hand is clearly shown 
in Fig.|^ (right panel): Z reflects the directionality in the polarization of this floating 
bond. Even if our results are still preliminary, they clearly show that the deviations of 
the effective charges from the crystalline value are noteworthy, and much larger than 
expected. 
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Table 1: Born effective charges (electronic) for some selected atoms in the a- Si super- 
cell. The effective charge tensor has been decomposed into an / = isotropic part 
(Zf~Tr(Z)/3), a Z = 1 antisymmetric part (Z^=(Z-Z'^)/2), and a / = 2 traceless 
symmetric part (Z'^=(Z-|-Z^)/2— Z^). The superscript T indicates transpose; atoms 
are labeled as in Fig. 1. 



4 Conclusions 

We have presented our results on the microscopic features of floating bonds in a-Si, 
which have been obtained with a maximally-localized Wannier functions approach. 
We confirm the conjecture that T5 defects are accompanied by well-defined delocal- 
ized states. Such states can be accurately characterized in terms of Wannier functions; 
a quantitative measure of delocalization is then provided by the the corresponding 
spreads. The delocalized states correspond to anomalous covalent bonds expanding 
over more than two atoms. The dielectric properties arc readily available as a byprod- 
uct of the Wannier analysis; we find strongly anisotropic effective charges that are 
significantly different from zero. 
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